Polymer-Enforced Crystallization of a Eutectic Binary Hard Sphere Mixture 
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We prepared a buoyancy matched binary mixture of polydisperse polystyrene microgel spheres of 
size ratio F = 0.785 and at a volume fraction of $ = 0.567 just below the kinetic glass transition. 
In hne with theoretical expectations, a eutectic phase behavior was observed, but only a minor 
fraction of the samples crystallized at all. By adding a short non-adsorbing polymer we enforce 
inter-species fractionation into coexisting pure component crystals, which in turn also shows signs 
of intra-species fractionation. We show that in formerly inaccessible regions of the phase diagram 
binary hard sphere physics is made observable using attractive hard spheres. 

PACS numbers: 82.70. Dd, 64.70.pv, 81.30.Bx, 61.66.-f 



Binary colloidal mixtures are valuable model systems 
for fundamental studies of crystallization [ - ']. Detailed 
predictions of their phase behavior exist for binary hard 
sphere (HS) systems of different size ratio T — Rs/Rl 
(where the Ri are the radii of small [S] and large [L) 
spheres, respectively). These show a sequence of spin- 
dle to azeotropic to eutectic phase diagram types with 
decreasing F, vanishing miscibility in the crystal phase 
for r < 0.85 and a huge variety of crystal structures for 
compounds [4-!)]. For zero miscibility eutectics the inter- 
species fractionation is expected to cause a pronounced 
slowing of nucleation [t^]. As the crystal phase becomes 
unstable for polydispersities above 6% [11, Vl] (10% in 2D 
[10]) an additional intra-species fractionation is expected 
for strongly polydisperse systems [ii;, 14]. Indications of 
polydispersity altered solidification dynamics have been 
observed in several experiments [15, 10] including sub- 
tle influences of the skewness of the size distribution [ ] 
but direct evidence of intra-species fractionation is still 
missing. While many of the predicted crystal structures 
have by now been observed [1, 2, 7, 18, 10] and also the 
general sequence of phase diagrams was confirmed [20], 
most parts of the phase diagram escape a detailed ex- 
perimental investigation. In particular, the exact loca- 
tions of phase boundaries in compound forming systems 
have not yet been determined. Crystallization kinetics 
has not been obtained and, except for a recent 2D-study 
[21], fractionation into coexisting L- and S'-crystals has 
not yet been observed. The main reason for this is the 
interference of crystallization with the glass transition at 
large volume fractions [22]. The second reason is the 
presence of gravity, which seemingly enhances the trend 
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to dynamically arrest the systems [ ] . As seen in other 
systems, gravity, moreover, may lead to differential sedi- 
mentation and inhomogeneities in composition [24, 25]. 

In the present paper we circumvent sedimentation by 
investigating a binary mixture of buoyancy matched mi- 
crogel particles. Exploiting observations made earlier on 
HS, copolymer micelles and binary HS [2(1-28], vitrifi- 
cation, on the other hand, is suppressed and crystalliza- 
tion is enforced by adding a short chained, non-adsorbing 
polymer. By adding polymer we move from HS to at- 
tractive HS (AHS), which, in principle, may significantly 
alter the phase behavior and crystallization kinetics [2!)- 

]. For the present case, however, we may state that the 
phase diagram type for binary AHS at this size ratio is 
identical to that of the corresponding HS system. For the 
first time we demonstrate a eutectic solidification process 
over the full range of compositions for a binary mixture of 
size ratio F — 0.785. We observe simultaneous precipita- 
tion of S- and i-crystals and give a preliminary account 
of their crystallization kinetics. As both species own a 
large polydispersity, we are, in addition, able to also ob- 
serve intra-species fractionation from the superposition 
of individual Bragg reflections. Further analysis suggests 
that the conversion efficiency is controlled by both the 
sample specific polydispersity and the polymer-induced 
mobility enhancement. 

Cross-linked polystyrene {PS) microgel spheres were 
synthesized by emulsion polymerization, cleaned, dried 
and re-suspended in 2-Ethylnaphtalene {2EN) as de- 
scribed previously [:!7]. In this good solvent the particles 
swell to hydrodynamic radii of Rh,s = (152 ± 4) nm and 
Rh.L — (193 ± 3) nm giving a size ratio of F;i = 0.785. 
The hardness of interaction was probed in oscillatory rhe- 
ological measurements following the procedure described 
in [-^S]. Assuming an inverse power potential U{r) oc l/r" 
we find n = 40±2 [ ] which is sufficiently steep to regard 
the particles as hard spheres. From form factor measure- 
ments the polydispersity index a — \J {R^) — (i?)^/ {R) 
is about 0.08 for both pure species. Sedimentation exper- 
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iments following [40] revealed the expected HS-like phase 
behavior. The freezing concentration was then identified 
with the ideal HS value and the melting points scaled 
accordingly to obtain $m,s = 0.528 and ^m,l = 0.531. 
Further, from the position of the first Bragg reflection of 
samples at coexistence the lattice spacing and from that 
the effective radii were obtained as Reff,L = (170.5 ± 2) 
nm and Reff^s = (126.2 ± 3) nm (effective size ra- 
tio ^eff = 0.74) [ ]. The volume fraction chosen for 
our experiments is just below the HS glass transition at 
= 0.573 ± 0.002 [42]. The phase diagram of the 
HS mixtures (i?l - R5) is shown in Fig. 1. To monitor 
the effects of attraction, a sample series {El - ElO) of 
eutectic composition {xs = 0.77) and increasing concen- 
trations, cp, of linear PS {Rg ~ 13.1 nm, — 133000 
g/mol; PSS GmbH, Germany) was prepared. Strongly 
attractive samples (^1 - Ab) of different molar fraction 
of S'-particles, xs, were prepared at cp = 12.54 g/1. Re- 
specting the volume occupied by the colloids, the polymer 
concentration in the free volume is cpjree = 43 g/1 [-. ], 
roughly 1.8 times larger than the overlap concentration 
c*p = 3Mw/i^TTNARl) = 23.5 g/1. 

Start of crystallization measurements (t — 0) was de- 
fined by taking a sample off the homogenizing tumbler 
and mounting it on a static light scattering goniometer 
(modified Sofica, SLS Systemtechnik, Germany, A = 405 
nm, 1 mW laser illumination, angular range 25° — 135°, 
resolution 1°). Structure factors, S{q,t), were obtained 
by time resolved measurements of the scattered inten- 
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FIG. 1: (Color online) Phase behavior of the pure binary 
mixture in the pressure-composition plane. The thin solid line 
gives the osmotic pressures for samples of $ = 0.567 ± 0.006 
and experimental size ratio — 0.785 as calculated using the 
fundamental measure theory equation of state of [ ] . Circles 
indicate the values for samples Rl - R5 with partial volume 
fractions cf>s = $s/$ of 0, 0.26, 0.57, 0.83 and 1 corresponding 
to molar fractions xs = Ps/{ps+Pl) of 0, 0.47, 0.77, 0.92 and 
1 (where pi = 'l>i/[(47r/3)-R| J is the particle number density). 
Thick solid lines are the phase boundaries adapted from [ ] 
for a HS mixture at Ft;, = 0.8. 



sity, I{q,t), for the concentrated dispersions, dividing 
by the appropriate form factors, Piq), obtained from 
a dilute sample and normalizing by the number den- 
sity ratio: S{q,t) = [Iconc{q,t)pdiiute]/[Piq)Pconc]- Here 
q ~ (47rnD/A)sin(6/2) is the scattering vector with 
npi denoting the suspension refractive index. Zero time 
structure factors, S{q, 0), only show the meta-stable melt, 
while at later times S{q, t) is a superposition of the con- 
tributions from fluid and crystalline phases. After two 
months most S{q, t) stopped evolving. The crystal scat- 
tering contribution, Sciq,t) was obtained following [4")]: 
Sc{q,t) = S{q,t) - l3S{q,0). The scaling factor /3 repre- 
sents the fraction of remaining melt and the total fraction 
of crystals Xc{t) = 1 — /3{t). In this approximation we 
neglect the small changes of melt and crystal composition 
and thus scattering power with time. 

The final S{q, oo) for the i?-samples in Fig. 2(a) are 
dominated by the scattering from the melt. R3 did not 
crystallize within two years and presented a single ex- 
ponential decay of its intermediate scattering function. 
All other samples showed isolated small crystallites. The 
final crystallinities Xc(oo) < 0.2 appear to be rather 
small. Presumably, values given by the lever rule apply- 
ing to the coexistence regions are further diminished due 
to polydispersity [13]. In each sample we could isolate 
individual narrow Bragg-peaks appearing at different q. 
Their average scattering vectors (dotted lines) reasonably 
well coincide with the melt maxima. The calculated os- 
motic pressures of the mixtures [ ] are below the eutec- 
tic pressure Pe = 28kBT{2RL)~^ at F = 0.8 [ ]. Hence, 
the observed phase behavior is consistent with theoreti- 
cal expectations but also with previous results on a 3% 
polydisperse mixture with F = 0.72 [ i]. 

The E samples in Fig. 2(6) show fluid order up to 
Cp — 9.05 g/1 and a minimum of the structural relax- 
ation time at cp = 8.12 g/1. For cp = 10.43 g/1 we 
observe crystallization, while for cp > 15.08 g/1 the sam- 
ples initially are amorphous again. Here the intermedi- 
ate scattering functions show indications of a plateau- like 
feature. E9 and ElO finally started crystallizing after two 
years. 

Fig. 2(c) shows the Sc{q, oo) of the A-samples with 
Cp = 10.43 g/1. Large numbers of small crystallites ap- 
peared within days and grew to final sizes ((0.4 — 0.8) 
mm) within several weeks. Neither additional peaks nor 
a shift in peak positions appears for the mixtures. This 
excludes the formation of compounds and of substitu- 
tional crystals. By comparison to the Sc{q, oo) of the 
two pure samples we identify A2 and ^3 to show crystals 
of both species, while AA contains 5'-crystals only. The si- 
multaneous occurrence of both crystal species shows that 
the systems are above the eutectic pressure (c.f. Fig. 1). 
In all mixtures the observed lattice constants equal those 
of the pure samples within experimental uncertainty. Us- 
ing the effective radii, i?e//.i, and the number densities 
inferred from the peak maxima we may calculate the vol- 



3 



10 20 30 40 50 




g/|jm' (//Mm"' g/pm' tid 

FIG. 2: (Color online) Scattering results. Curves are shifted for clarity, (a) Final structure factors of _R-samples at different 
compositions as indicated. $(_R1) — 0.54, <1>(_R2, i?3, i?4) = 0.567, "I>(i?5) = 0.505. Dashed lines indicate location of Miller- 
indexed reflections averaged over several individual crystallites. (6) Final structure factors for the iJ-samples of xs = 0.77 = xe 
and $ = 0.567 with indicated amounts of added polymer, (c) Crystal structure factors of A-samples with location of Miller- 
indexed reflections of the pure samples averaged over several individual crystallites (dashed lines), (d) Temporal development 
of the fraction of crystals for the three mixtures A4 - A2. Filled symbols denote 5-crystals, open symbols denote L-crystals. 



ume fractions for both coexisting crystal species. In all 
samples we obtained values of = 0.73±0.02. This is re- 
markably close to the maximum packing fraction obtain- 
able in a close packed structure. Similar compressions so 
far have only been reached at this size ratio between col- 
loids and depletant, when additional gravitational com- 
pression of sedimented crystals occurred [ '"]. 

From time resolved measurements the fraction of crys- 
talline material, X{t), is obtained by integrating the in- 
tensity in the corresponding peak regions: XL,s{i) — 
c J dqSc{q), choosing c such that ^^(t) + Xs{t) ~ 
Xc{t) = 1 - I3{t). Fig. 2(d) shows that at xs = 0.92 
about 20% of the 5-particles form S'-crystals already at 
early times which grow for about a month; at xs — 0.77, 
5% of each species form crystals simultaneously; the final 
crystallinity is reached after 10 days; for xs — 0-47 start- 
ing again from about 5% of each species, S'-crystal con- 
version accelerates after L-crystals have nearly reached 
their final crystallinity. In all samples the degree of crys- 
tallization Xc{co) ~ 0.4 — 0.6 is considerably enhanced 
over that reached in the i?-samples. Moreover, at this cp 
the final Xs of 0.42, 0.4, 0.3 and 0.2 roughly scale with 
the number fractions of 1.0, 0.92, 0.77 and 0.47, show- 
ing that about 40% of the S'-particles are able to solidify 
via in-species fractionation. Similarly, for the L-particles 
we find Xl to be 0.13, 0.34 and 0.6 for (1 — p) being 
0.23, 0.53 and 1.0, i.e. more than half of the particles 
crystallize. 

Summarizing, we find that the addition of non adsorb- 
ing polymer to a buoyant binary mixture allows investi- 
gations in regions of the phase diagram, in which under 



gravity much more complicated scenarios occur [2-5, 32] 
and which in pure HS usually are blocked by a kinetic 
glass transition. For the particular system investigated 
here the reference system without polymer shows fluid- 
solid coexistence or pure fluid. The observed phase be- 
havior is consistent with the theoretically expected eutec- 
tic [4-S]. By contrast, [L + S) coexistence is obtained, 
when the pressure, P, surmounts the eutectic pressure 
Pe{xetCp_e), corresponding to the eutectic composition 
and a polymer concentration 9.05 g/1 < cp^E < 10.43 g/1. 
Recent work on one component HS [35, 36] showed that 
both and Pp shift to lower values upon increasing 

the attraction. The present work suggests that for bi- 
nary AHS Pp drops faster than P(<I> = 0.567). However, 
theoretical predictions for composition dependent equa- 
tions of state and the phase diagram of binary attractive 
HS are eagerly awaited to quantitatively rationalize our 
observations in a plot analogous to Fig. 1. 

Being above Pe, our mixtures should convert com- 
pletely to S- or L-crystals. Instead, we find a pronounced 
amorphous background signal and broad pyramid-shaped 
peaks composed from a superposition of a large number 
of individual Bragg reflections located at slightly differ- 
ent q. This supports theoretical expectations that con- 
version proceeds via combined inter- [8] and intra-species 
fractionation [J ]. In both cases the conversion efficiency 
should be correlated with the possibility to form com- 
bined density and composition fluctuations. The addi- 
tion of polymer generally enhances the dynamics of such 
fluctuations [7] as long as the system remains untrapped 
by the attractive glass. However, conversion efficiency 
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also relies on the availability of similar sized partner par- 
ticles to form a stable fluctuation. We here observe an 
absence of L-crystals at large Xg (already seen in [19, 25]), 
an increased initial formation of iS-crystals in A4, their 
delayed formation in A2, and an overall incomplete con- 
version. This suggest that depletion of the minority com- 
ponents (initial and/or in time) and their enrichment in 
the melt (by formation of majority crystals) govern their 
fluctuation spectrum and conversion efficiency. Particles 
from the wings of the size distributions are simply too 
diluted to crystallize. 

In this study we have recovered important aspects of 
HS-crystallization in an AHS system. In future the con- 
trolled use of depletion attraction in formerly inaccessi- 
ble regions of the phase diagram may allow detailed and 
quantitative studies of the phase behavior at other size 
ratios, kinetics of compound formation or other interest- 
ing problems of condensed matter physics. 
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